Dirac Nodes and Quantized Thermal Hall Effect in the Mixed State of d-wave 

Superconductors 
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We consider the vortex state of d-wave superconductors in the clean hmit. Within the linearized 
approximation the quasiparticle bands obtained are found to posess Dirac cone dispersion (band 
touchings) at special points in the Brillouin zone. They are protected by a symmetry of the linearized 
Hamiltonian that we call Toirac- Moreover, for vortex lattices that posess inversion symmetry, it is 
shown that there is always a Dirac cone centered at zero energy within the linearized theory. On 
going beyond the linearized approximation and including the effect of the smaller curvature terms 
(that break Toirac), the Dirac cone dispersions are found to acquire small gaps (« 0.5KTesla~^ 
in YBCO) that scale linearly with the applied magnetic field. When the chemical potential for 
quasiparticles lies within the gap, quantization of the thermal-Ifall conductivity is expected at low 

temperatures i.e. — n{ ) with the integer n taking on values n — ±2,0. This quantization 
could be seen in low temperature thermal transport measurements of clean d-wave superconductors 
with good vortex lattices. 



I. INTRODUCTION 

Since the experimental verification of the dj.2_y2 na- 
ture of superconductivity in the cuprate materials |^] 
there has been much activity in studying the physics of 
quasiparticles in a d-wave superconductor. In contrast 
to the s-wave case, the d-wave superconducting gap van- 
ishes at points on the Fermi surface leading to low energy 
quasiparticles that behave like massless Dirac fermions. 
One question that is of great interest is the behaviour of 
these quasiparticles in the mixed state of superconduc- 
tors - a problem that is both theoretically rich as well as 
of relevance to explaining several different experiments. 
When the vortices are arranged in a lattice, the situa- 
tion is reminescent of the Hofstadter problem of charged 
particles in a magnetic field, subject to a periodic poten- 
tial, where much interesting physics is known to emerge. 
Here, however, the periodic lattice is generated by the 
magnetic field, and, as a result of flux quantization, is 
always in commensuration with the field. The quasipar- 
ticles, which do not carry a definite charge, also couple 
differently to the field - in fact they only couple to the 
superflow arising from the combined effect of the mag- 
netic field as well as the vortices. Furthermore, there is a 
statistical interaction between hc/2e vortices and quasi- 
particles, which causes them to change sign on circling 
a vortex. As we shall see, these properties will intro- 
duce interesting new features into the problem. Experi- 
mentally, there have been several probes of quasiparticle 
behaviour in a magnetic field, including measurement of 
thermal hall conductivity Kxy and low temperature 
longitudinal thermal conductivity Kxx and specific 
heat 1^ . A satisfactory explanation of the results of these 
experiments in terms of da,2_y2 quasiparticles would be 
strong support for the point of view that the supercon- 
ducting phase of the cuprates is a conventional d-wave 
superconductor. On the other hand, if such an explana- 



tion proves elusive, it may well point to some additional, 
and perhaps exotic, physics even in the superconducting 
state of the cuprates. 

This problem of Dirac quasiparticles in the mixed state 
has been considered by several authors. Gorkov and 
Schriffer |Q, and Anderson |^ proposed a Landau level 
like spectrum. In the latter work, this was derived by 
mapping the problem onto a Dirac particle in a uniform 
magnetic field, under the assumption that the superfiow 
may be neglected. In Morita et al. studied a lat- 
tice model of the problem, first in the same approxima- 
tion as 1^, and then numerically with the superfiow in- 
cluded. Topological aspects of relevance to this problem 
were pointed out. 

A different approach was taken by Franz and Tesanovic 
who considered the problem within the linearized approx- 
imation and introduced a ficticious U(l) gauge field to 
implement the statistical interaction between the vortices 
and quasiparticles (Franz- Tesanovic transformation) [p^ . 
This led to a numerical evaluation of the resulting quasi- 
particle band structure, which was extended by the de- 
tailed study of Marinelli et al. , Vafek et al. and 
Kallin et al. H. 

Here too we begin by considering the the problem of d- 
wave quasiparticles in the vortex lattice state within the 
linearized approximation. In contrast to some earlier ap- 
proaches, we focus on identifying the symmetries of the 
problem and their consequences for the spectrum of the 
linearized Hamiltonian. Our results are easily stated - the 
Hamiltonian, as a consequence of linearization, posesses 
an additional symmetry (which we call Tuirac) that pre- 
serves Dirac cones at certain special points in the Bril- 
loin zone. The energy dispersion, in the vicinity of these 
points, is that of a massless Dirac particle. If the vortex 
lattice posesses inversion symmetry, then there is a Dirac 
cone centered at zero energy. 

While these results are in general agreement with the 
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results obtained in earlier numerical work [|l0|-|l3[ we 
point out out a subtle feature of the Franz-Tesanovic 
transformation in the linearized approximation, which 
could yield spurious gaps in numerical simulations and 
which we believe to be the cause for discrepancies from 
the results derived here. 

We then proceed to consider the effect of the curvature 
terms, that were dropped when linearizing the Hamilto- 
nian. These terms arise, for example, from the parabolic 
nature of the electron dispersion. In the parameter range 
of interest they may be considered as small perturbations 
on the linearized Hamiltonian. However, as pointed out 
in [0. they are crucial to generating a non- vanishing 
thermal Hall response. Since the curvature terms break 
the symmetry Toirac that protects the Dirac nodes, they 
give rise to a small gap at these nodes, and the energy 
dispersion in their vicinty is now that of a massive Dirac 
particle. Thus for temperatures that are smaller than 
this gap scale, the curvature terms can have a qualita- 
tive effect on the properties of the system. 

In addition to information about the quasiparticle 
spectrum, it is in some cases possible within our ap- 
proach to derive consequences for low temperature ther- 
mal transport. To this end, we first recall that mas- 
sive Dirac particles in two dimensions exhibit a quan- 
tized Hall effect when the chemical potential lies in the 
gap In the presence of a vortex lattice with 

inversion symmetry, the chemical potential for the su- 
perconductor quasiparticles will lie at the center of the 
gap induced by the curvature terms, if we ignore Zeeman 
splitting. As a result, a quantized thermal Hall conduc- 
tance of Kxy/T = ±1/2 (in appropriate units, at low 
temperatures) is expected from each of the four nodes, 
and can give rise to two scenarios Kxy/T = ±2 or 0, while 
Kxx/T ^ in both cases. First, if the contribution of all 
four nodes is of the same sign, Kxy/T = ±2, we have 
a nontrivial quantized thermal Hall conductance that is 
expected to be attained for temperatures smaller than 
the gap. This situation is topologically identical to (i.e. 
has the same edge content as) a pure dx2-y2 + i dxy su- 
perconductor, that is also known to exhibit a quantized 
thermal hall effect |l|]. A magnetic induction of 
such a pairing symmetry in the cuprates was proposed 
in jl^ 1^ ; here we will provide a concrete realization of 
these general ideas, and layout the route to calculating 
physical parameters such as, for example, the size of the 
energy gap. The second scenario is when the Hall con- 
ductances from the different nodes cancel leading to a 
Kxy/T = 0, that is topologically identical to a dx2-y2+is 
superconductor, or any other thermal insulator. Which 
of these scenarios is realized is a function of vortex lattice 
geometry and the anisotropy of the Dirac dispersion in 
the homogenous superconductor. In order for the quanti- 
zation to be visible experimentally the energy gap needs 
to be larger than the Zeeman splitting, which plays the 
role of chemical potential for the quasiparticles . Since 



the energy gap is also found to scale linearly with mag- 
netic field, and is roughly of the same magnitude as the 
Zeeman energy, the question of which one is larger in any 
particular material is a detailed quantitative issue. Here 
we will demonstrate how fairly simple numerical calcula- 
tions within our theory can predict in a given situation, 
which of the scenarios described is realized, as well as 
provide a quantitative estimate of the energy gap and its 
dependence on various physical parameters. 

The rest of this paper is organized as follows. In Sec- 
tion II, we begin by laying out our assumptions and 
then deriving the Bogoliubov-de Gennes equation that 
describes d-wave quasiparticles in the mixed state. The 
quasiparticles are shown to couple to the superflow gen- 
erated by the combined effect of the vortices and the 
magnetic field. In addition, they acquire a Berry phase 
of (-1) on circling a hc/2e vortex. We then derive the 
linearized approximation that describes the low energy 
quasiparticle excitations in the magnetic field range of 
interest. These linearized equations describing d-wave 
quasiparticles in a vortex lattice are analyzed in Section 
HI. For purposes of clarity, we first consider a vortex lat- 
tice of hc/e (double) vortices. The simplifying feature 
here is that the Berry phase terms are absent and we 
only have to contend with the effects of the superflow. 
A symmetry of the linearized Hamiltonian that protects 
the Dirac nodes is identified, and the role of inversion 
symmetry in maintaining a Dirac node at zero energy is 
described. We then tackle the physically more interesting 
case of a vortex lattice of hc/2e vortices. Here, a similar 
though more involved analysis obtains for us the same re- 
sults as for a vortex lattice of hc/e vortices. A check on 
these arguments within perturbation theory is detailed 
in Appendix A. In Appendix B, we turn to some sub- 
tle aspects of the Franz-Tesanovic (FT) transformation 
for the linearized problem, especially the issue of gauge 
invariance under different FT transformations, and the 
need to properly regularize the linearized theory in order 
to obtain a faithful representation of the problem. An al- 
ternate argument for the existence of Dirac nodes in the 
linearized problem for the hc/2e vortex lattice case, that 
docs not utilize the Franz-Tesanovic transformation, is 
also presented in this appendix. A detailed comparison 
of our results for the linearized theory against earlier nu- 
merical work is presented in Appendix C. In Section IV, 
we go beyond the linearized approximation by including 
the effect of curvature terms and discuss implications for 
low temperature heat transport. A numerical calculation 
of the gaps induced by the curvature terms for the hc/e 
square vortex lattice case is also presented. Finally, in 
Section V we conclude with some brief comments on the 
effect of vortex disorder and a comparison of our theoret- 
ical expectations against available experimental data. 
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II. THE BDG EQUATIONS FOR D-WAVE 
QUASIPARTICLES IN THE MIXED STATE 

In this section we derive the equations governing quasi- 
particles in a dx^-y^ superconductor, in the presence of 
a vortex lattice. We begin by detailing the assumptions 
and approximations that we will make in what follows. 



A. Assumptions and Approximations 

(a) Existence of quasiparticles: The systems that we 
will mainly be interested in are the cuprate high tempera- 
ture superconductors, that are known to have dx2-y2 gap 
symmetry. We assume that this superconducting state is 
otherwise conventional, in particular that there are well 
defined quasiparticle excitations in this phase, for which 
there is experimental support from angle resolved pho- 
toemission studies . For inhomogenous situations, such 
as the vortex lattice state, we assume that the quasi- 
particles are governed by an appropriate Bogoliubov-de 
Gennes type equation. 

(b) Neglecting Vortex Core Contributions: The vor- 
tex core is the region around the center of the vortex 
of size ^, the coherence length, where the magnitude of 
the order parameter is significantly supressed from its 
bulk value. At fields much smaller than iJc2, the vortex 
cores of extreme type II superconductors are significantly 
smaller than the separation between vortices. This is the 
situation, for example, in optimally doped YBa2Cu306.9 
(YBCO) over the accessible field range. The vortex cores 
have a size ^ 15 A while the inter- vortex separation is 
of order 500 A in a 1 Tesla field. Since the vortex cores 
take up so little of the sample area (0.1% in this exam- 
ple) we neglect the modulation of the order parameter 
magnitude, while retaining its phase variation. 

(c) Perfect Vortex Lattice: We will make the assump- 
tion that the superconductor is clean and the vortices 
are arranged in a perfect lattice. Of course, in any real 
situation, vortex disorder is expected to be present. In 
some cases it may even be so large as to destroy the long 
range positional order of the lattice, in which case, of 
course, the perfect lattice approximation will not be a 
good starting point - for example, in BSCCO in mag- 
netic fields of around one Tesla, neutron scattering stud- 
ies indicate that an ordered vortex lattice is absent . 
However in YBCO, Bragg spots from the vortex lattice 
are seen for which the starting point of a perfect vor- 
tex lattice may be more justified. Despite the existence 
of vortex disorder, in what follows we shall proceed with 
the assumption of a perfect vortex lattice as it provides 
us with a theoretically tractable starting point. Besides, 
there is some evidence from experiments and theory 
p3[ that the scattering from vortex disorder at low tem- 
peratures is rather small, and neglecting its effect may be 



permissible at a first approximation. Extending the the- 
ory to include the effects of weak vortex disorder is left 
for future investigation, although we briefly return to the 
topic of vortex disorder in Section V, while considering 
the stability of our results. 



B. Bogoliubov-de Gennes equations for d-wave 
superconductors in the mixed state 

Consider a two dimensional d-wave superconductor in 
the pure state, which may be described by the model 
Hamiltonian 

T-(-pure = I ^ct^e(p)cp^ 4- i(A(p)c];^T^^-clp^, +h.c.) 

(1) 

where c^^ is the creation operator for an electron of 
momentum p and spin projection a and t^^i are elements 
of the unit antisymmetric matrix: 



1 

-1 



(2) 



Here we consider a parabolic dispersion for the electrons: 

1 



2m 



P 



Ep 



(3) 



and a d^y symmetry of the gap function (equivalent to 
da;2_j,2 rotated by 45^^) which could be taken to have, 
near the Fermi momentum, the functional form: 



A(p) 



Ao 



(4) 



. We now consider the effect of a magnetic field applied 
along the z axis, that gives rise to vortices. It is conve- 
nient at this stage to define the 'd' operators as: 

dip — Cp'^ 



d2p = c 



'Pi 



here the spin projection is taken in the direction of the 
magnetic field. The Hamiltonian for a layer of supercon- 
ductor in the mixed state can them be written as: 



n = 



(fr 



(4(r) d\{r))HBaG[^l^l^h+nz^ 



(5) 



where 

HbcLG 



e{-iV-eA{r)) e'^ A{-iV)e-'^ 
e-*^ A(-iV)e''^ -e(iV-eA(r)) 



where A is the vector potential V x A{r) = 
(j){r) is the phase of the order parameter; and 



(6) 

B{r)z and 



3 



n 



Zeeman 



d'^rB{f)[4{r)di{r) + 4(^)^2(^)1 



(7) 

. The Hamiltonian written in terms of the 'd' particles 
contains no anomalous terms which implies that the to- 
tal number of 'd' particles is conserved. This is simply a 
reflection of the fact that the z component of the spin is 
conserved in this situation. Thus, the density of 'd' par- 
ticles, poir) = {d\{'r)di{r) + dl{r)d2{'r)) is proportional 
to the spin density of the quasiparticles along the z axis, 
and so appears in the expression for the Zeeman coupling 
to the magnetic field. The uniform part of the magnetic 
field density will therefore behave as a chemical potential 
for the 'd' particles. In what follows, we shall not retain 
the Zeeman term, but rather comment on its effect at the 
very end. 

In equation (|^) the phase variation of the order pa- 
rameter, induced by the vortices, has been written in the 
form e'a' A(— zV)e~*^ . This preserves gauge invariance 
and is consistent with the formulae in |^,^. The vari- 
ation of the magnitude of the order parameter has been 
dropped, as discussed earlier. We shall take as given the 
magnetic field distribution B{f) and the phase variation, 
which must satisfy the equation 



V X (V0) = 27r ^ (5(2) (f-i?,)^ 



(8) 



where RiS denote the positions of the vortices and the 
sum runs over all vortices in the lattice. 

Due to the conservation of the 'd' particles, and the 
fact that they do not interact with each other at this 
level, we can write a wave equation which contains all the 
physics of this many body system. This wave equation 
is the Bogoliubov-deGennes equation: 



H 



BdG 



u{r) 
v(r) 



= E 



u{r) 
v{r) 



(9) 



where HbcLG has been defined in equation (^ and can 
be written in the compact form: 



HBdG = e(-iV-eA(r)cr,)cr,+{e*— A(-iV)e 



2 <T 



u 



The transformed Hamiltonian takes the simple form, 



H'saa = UHBdoU-^ = e(-iV + P.(r)^,)o-, 



A(-iV)cr^ 

(11) 



where Pg = ^Vt/i — eA a, gauge invariant quantity, is 
the mechanical momentum carried by each member of 
Cooper pair at point r. We will sometimes refer to this 
quantity as the superflow, though this terminology is not 
quite accurate. Given this definition of Ps, it is easily 
seen that 



V X A- = [{i/,^5(f - - eB{r)V- 



(12) 



In the presence of elementary ^ vortices, it must be 
noted that the unitary transformation ( |lO| ) is not single 
valued, but changes sign on circling an odd number of 
such vortices, since it depends on the half angle 0(r)/2 
that winds by an odd multiple of tt. Thus the trans- 
formed quasiparticle wavefunctions, {u' v')'^ defined by, 



u'ir) 
v'{r) 



= U 



u(r) 
v{r) 



(13) 



are not single valued, but change sign on circling an 
odd number of ^ vortices. Methods for handling this 
statistical interaction between quasiparticles and vor- 
tices, are described in the next subsection. 

Thus, writing the Bogoliubov-deGennes equation in 
this gauge clarifies how the quasiparticles interact with 
the magnetic field - which is as summarised below. 

• The quasiparticles couple to the superflow set up by 
the combined effect of the vortices and the magnetic 
field, namely to the combination Pg = — eA, 
and not individually to the phase variation or the 
vector potential. 

• Pick up a (-1) Berry phase factor on circling a ^ 
vortex. 

• Also interact with amplitude variations of the order 
parameter, which we have dropped. 



-h.c.} 



The crs are the 2x2 Pauli matrices in the usual repre- 
sentation. (In this language, the Zeeman term takes the 

form H Zeeman = ^-B(r)l). 

It is convenient to make a gauge transformation to 
eliminate the phase variation from the order parameter 
(London gauge) which may be affected by the unitary 
transformation: 



(10) 
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FIG. 1. The gap structure for a d^y superconductor with 
nodes along the coordinate axes. 

C. The Linearized Approximation 

In the limit of low temperatures, and magnetic fields 
much smaller than Hc2, we can restrict our attention 
to the low energy quasiparticle excitations near the gap 
nodes of the d-wave superconductor. We first recall the 
linearized approximation for the pure case, which gives 
us the Dirac dispersion of the nodal quasipartcles. The 
Hamiltonian for quasiparticles with momentum 'p' is 
Hp = e(j))<T z + A(j))(Tx ■ For low energy excitations in the 
vicinity of node 1 (see Figure this can be expressed in 
terms of p — ppx + Sp (where Sp is measured from node 
1), as 




(14) 

where vp — pp/m, wa — Aq/pf- The linearized ap- 
proximation consists of dropping the terms within brack- 
ets, in which case we obtain an anisotropic Dirac disper- 
sion for the quasiparticles near the nodes:- 

El{p) = vp^iSp^y +a~^{Spy)^ 

where a — vp/vA, is the anisotropy of the Dirac disper- 
sion. Following [0 we estimate the temperature scale up 
to which this approximation can be trusted. We require 
that the value of the quadratic terms at the typical mo- 
menta of the thermally excited quasipartices be smaller 
than the linear contribution. For YBCO, this tempera- 
ture turns out to be Tj^ax ~200Kelvin p^ . 

We now consider a similar low energy approximation 
for the vortex state of the dxy superconductor. First, we 
describe a proceedure to handle the non-single valued na- 
ture of the wave functions tp' = ^ ^ '^^ equation ( [T^ ) 

which is called the Franz- Tesanovic transforma- 
tion, although we adopt a slightly different approach to 
its derivation here. 

Franz- Tesanovic Transformation: Recall that we want 
to solve the eigenvalue equation; 

H'Bdoi'' = Ei;' (15) 

with the condition that tp'{x,y) is not single valued, but 
acquires a negative sign on circling an odd number of ^ 
vortices. To implement this condition, we write ip' as a 
product of a fixed function that is multiple valued which 
precisely builds in the required sign changes (<I>(a;, ?/)), 
and a single valued wave function tpFT{x,y). Thus, 

ip'{x, y) = $(x, y)ipFT{x, y) (16) 



, and the preceeding eigenvalue problem can be for- 
mulated as an equivalent problem for the single valued 
wavefunction ^pt- In order that this is also a hermitian 
eigenvalue problem, we choose $(a;, y) from the class of 
functions: 

where the qi are odd integers, and we have used the com- 
plex coordinates z = x + iy, z = x — iy. The product 
runs over all || vortices, located at Zi = Xi + iyi. Note 
that this function is pure phase i.e. |${q}.(a;, y)| = 1 [ p7[ . 
The choice of the odd integers qi is arbitrary and repre- 
sents a gauge degree of freedom. Clearly, physical results 
cannot depend on this choice, though in practice we may 
work with a particular set of {qi\ that is convenient for 
calculation. Due to its singular nature, this transforma- 
tion needs to be handled with some care especially in 
the linearized theory, a point which we will return to in 
Appendix B. 

Inserting ^ in (|l|), we have: 

HpTlpFT — EtppT (18) 

HpT ^ e{p + d+ Ps(Tz)(Tz + /:^{p + a)(T^ (19) 

where a(a;, y) is a real vector field given by: 

a = -iW log i>{,} (f) (20) 
h r — ri 

2 ^-^ \r — n r 

i ' ' 

which implies 

V X a{r) = fnT'^qi6^'^\f~ fi) 

i 

that solenoids of flux nqi of the a gauge field have been 
attached to the vortices. The sign change of the quasipar- 
ticles on circling a unit vortex is now accounted for by the 
Aharonov-Bohm effect arising from this solenoid of flux. 
Thus, a flcticious U(l) gauge fleld has been invoked to 
handle the (-1) phase factors acquired by a quasiparticle 
on circling a vortex. Clearly, this is a highly redundant 
description - in principle the U(l) gauge field a can ac- 
count for any phase factor, which is here being restricted 
to just ±1. The minimal choice that would take care of 
just these two factors, is an Ising (Z2) gauge field, which 
however requires a real space lattice for its formulation 

0- . n 

Linearization for the Vortex State: Following |14[ , we 
consider low energy excitations near the nodal points, 
and neglect the effect of inter-node scattering. Then, we 
can expand the wavefunction as: 

^^y(f) = e^'^^^ij^ir) + e-^'=^"i^i(f) 

+e'''^yMr) + e-'^'^yij^ir) (22) 
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where the functions tpi (i=l,T,2,2) are considered to be 
slowly varying on the scale of kp^. The problem then 
reduces to solving, at each node: 



(23) 



where the Hf^ represent the linearized part of the 
Hamiltonian dlq): 



FT 



(24) 

= vf{{Px + ax)(Tz + a^^ipy + ay)(Tx + x ■ Ps{r)} 



i?2 - -H^ 



FT 



(25) 

= VF{{Py + ay)(Tz + a~'^{px + ax)(Tx + V ■ Ps{r)} 

with p = — iV. The remaining part of the Hamiltonian, 
Ai7, arises from the curvature of the electon dispersion 
and of the gap function and is given by: 



2m 



p^^m^^ + ^ip 



2 m 



-■^{Pi +a.i,P2 + a2}o-x 
^Pf 



(26) 



where {,} denotes antisymmetrization; {A, i?} — AB + 
BA. 

It may now be argued that the linear piece Hf^, dom- 
inates over the curvature terms AH for magnetic fields 
much smaller than Hc2- Let us denote the typical sepa- 
ration between vortices by d oc 1/^/B. While the linear 
part of the Hamiltonian has terms of order Ei — hvp/d 
(50K Tesla-5 in YBCO) or of order E2 ^ hvA/d (3K 
Tesla~2 in YBCO), the largest terms in AH are of order 
Ef/Ep (0.5K Tesla~^ in YBCO), so long as we assume 
that the quasiparticle trajectories do not come too close 
to the vortex core. Thus, in the magnetic field range of 
interest of a few Tesla, we can at a first approximation 
neglect the curvature terms AH and simply solve the lin- 
earized problem. Subsequently we will include the effect 
of these curvature terms (they are crucial to producing a 
finite thermal Hall signal), and verify that their effect is 
indeed a perturbation on the linearized problem. 

We are therefore led to consider the linearized problem, 
say at node 1: 

H[^i} = Ei; 

H['^ = vpiipx + ax)(Tz + {Py + ay)(Tx + Ps^(f)l} (27) 

We have already noted, that in addition to the interac- 
tion with the amplitude variation of the gap (which has 
been neglected) the quasiparticles interact with the vor- 
tices via the superflow Pg and the Berry phase of (-1) 
acquired on circling a fundamental vortex. For the sake 
of clarity, it is helpful to first consider a situation where 
the Berry phases are inactive, and the quasiparticles only 
see the superflow. This is formally accomplished by con- 
sidering the case of a vortex lattice of (double) vor- 
tices. Then, there are no nontrivial Berry phase factors 



and hence the a fields (which were inducted to keep track 
of these phase factors) , may be dropped. Equivalently, it 
may be noted that for the case of vortices, the trans- 
formation ( p^ ) is single valued. Below, we will analyze 
the linearized problem in the presence of a ^ vortex lat- 
tice. Armed with this understanding, we then consider 
the physically relevant case of a lattice of || vortices. 
Although our conclusions from the — vortex lattice case 
go over unchanged, the reasoning there is a little more 
involved. 



D. Vortex Lattice of — vortices in the linearized 

e 

approximation: 

Consider the linearized hamiltonian at node 1, for this 
conceptually simpler case of a vortex lattice of ^ vor- 
tices. 



Hi = VF{PxCrz+a ^pyO-^ + Psx{r)l} 



(28) 



where the superflow Pg = {Psx , Psy) is given by the gauge 
invariant combination, 

Ps = [iv0 - eA] 

and satisfies, 

V xP,^[{hY,5''{r-Ri))-eB{r)]z (29) 

i 

where the sum runs over all the vortices in the lattice, at 
positions Ri. This differs from equation (|l2|) in that the 
strength of the vortices is now doubled. For a periodic 
lattice of vortices, Pj, being a physical quantity is also 
a periodic function with the same period as the vortex 
lattice. Formally, this can be seen by noting that:- 



^1 



/unit cell 
from flux quantization. Thus, 



(fr{V X Ps) =0 



(30) 



Ps-dl = 



on the boundary B of the unit cell, which is consistent 
with a periodic superflow Pg. 

Consequently, the Hamiltonian Hi (|2^ ) is that of a 
Dirac particle in the presence of a periodic scalar poten- 
tial (played by Psxir)). Cearly, a band structure will 
result, and the eigenstates can be labelled by the band 
index, and the crystal momentum k that takes values 
within the Brillouin zone. What follows is a symmetry 
analysis of the spectrum of Hi . One of the main issues 
we address is whether the Dirac node at zero energy, 
that is present for the free Dirac Hamiltonian (and cor- 
responds to the node of a dxy superconductor in the pure 
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state) survives in the presence of the periodic superflow 
term Psxir)- Our results are as follows. Dirac nodes 
(band touchings) survive in the spectrum of Hi at those 
points in the Brilloin zone that are invariant under the 
k ^ —k transformation. The symmetry that protects 
these nodes, which we call Tjjirac, is obtained as a conse- 
quence of the linearization. Further, we find that there is 
a Dirac node centered at zero energy if the vortex lattice 
posesses inversion symmetry. 



III. DIRAC NODES IN THE LINEARIZED 
PROBLEM: 

We begin by analyzing the conceptually simpler case of 
a vortex lattice of hc/e vortices, where the Berry phase 
factors are absent, and we have only the coupling of the 
quasiparticles to the superflow to deal with. We then 
tackle the hc/2e vortex lattice case along parallel lines. 

A. Vortex Lattice of — Vortices 

e 

We will consider the linearized Hamiltonian for node 1 
Hi = vfIPxCTz + a'^pycr^ + Psxirjl} 



which is equation (^8|) of the previous section, where 
Ps(r) for an arbitrary lattice of vortices, is a periodic 
function with the same period as the lattice and satisfies 
equation (|2^) . This gives rise to a band structure for the 
quasiparticles, and the eigenstates are labelled by a band 
index and a crystal momentum [k), which takes values 
within the Brillouin zone. We consider for generality an 
oblique vortex lattice, which posess no additional sym- 
metries. The Brillouin zone for such a vortex lattice is 
depicted in Figure ||. 

A Symmetry of the Linearized Hamiltonian: Consider 
an eigenstate ipj: of the linearized Hamiltonian Hi with 

eigenvalue E and crystal momentum k: 



Then the transformed wavefunction (f'_j: = Tip'^ is also 

an eigenstate with energy E but crystal momentum — fc, 
where r is the dimension two antisymmetric matrix: 



T = 



1 

-1 



This can be verified by noting that the linearized Hamil- 
tonian is taken to itself under the transformation: 



- tHIt = Hi 
where we have used the identities 



(31) 



= -1 
TCr;*X = (Ti 

Since this transformation is formally equivalent to the 
time reversal operation for Dirac particles, we will call it 
TDirac, although, as explained below, it is distinct from 
the physical time reversal transformation for this prob- 
lem. This symmetry ensures that states with crystal mo- 
mentum k and —k have the same energy. For those points 
in the Brillouin Zone that are taken to themselves under 
time reversal (modulo a reciprocal lattice vector: k = —k 
[mod G)]) there is a degenerate pair of states t/jj; and tV"^- 
These states are orthogonal, from the antisymmetry of 
T. Hence we find a degenerate doublet at these special 
points in the Brillouin Zone, which are shown in Figure 
g as the r (fc = 0), A, B and M points. The spectrum 
at these special points of the Brillouin Zone is composed 
entirely of degenerate pairs, and as we will see shortly, 
on moving a little bit away from these points in crystal 
momentum, the states split and give rise to the energy 
dispersion of a massless Dirac particle. 




• r 








FIG. 2. The Brillouin zone for an oblique lattice of vortices. 
The points of the Brillouin zone that are taken to themselves 
under time reversal symmetry are shown as F, A,B and M 
points. At these points degenerate doublets and hence Dirac 
cones are expected within the linearised theory for all lattices 
of — vortices. The same holds for lattices of vortices. 

e 2e 

The symmetry operation Tuirac that operates on the 
quasiparticle excitations at a single node is distinct from 
the physical time reversal operation that would trans- 
form states at one node into states at the opposite node. 
Rather, Toirac is obtained as a symmetry as a conse- 
quence of linearizing the electron dispersion, and it is 
easily seen that the curvature terms, such as for example 
^p^eXz, violate this symmetry. Thus, although T^nrac is 
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not a symmetry of the entire problem, since the curvature 
terms that violate it are so small, it is still a good ap- 
proximate symmetry. For the linearized theory of course, 
it is an exact symmetry. 

Dime Cones from Degenerate Doublets: As we move 
away from the special points in the Brillouin zone at 
which degenerate doublets are found, the crystal momen- 
tum splits these states and gives rise to a Dirac cone 
dispersion. This is most easily appreciated by analogy 
with the pure case. There, the free Dirac Hamiltonian, 
Hfree = vpPxO'z+VAPyO'x, has a degenerate pair of states 
at the center of the Dirac cone (at p = 0) given by the 
two constant spinors that are eigenstates with zero en- 
ergy. Moving away in momentum, which is exactly anal- 
ogous to turning up the Zeeman field on a spin one-half 
particle, the states are split into ±E pairs, with the split- 
ting being linearly proportional to the momentum. This 
give rise to the massless Dirac dispersion that we expect 
for this Hamiltonian. In a similar way, we find here that 
Dirac cones arise centered at the special points at which 
the degenerate doublets are present. 

Consider a pair of degenerate states, ^^{r) and Ti/j*{r), 
at one of these special points in the Brillouin zone. The 
effect of moving away from this point by crystal momen- 
tum Sk can be accounted for by adding the piece: 



5Hsk = vpSkxCTz + v^SkyCTx 



(32) 



to the Hamiltonian and leaving unchanged the boundary 
condition for the wavefunction. If the deviation in crystal 
momentum is small (compared to the reciprocal lattice 
vectors), this additional piece will barely mix the different 
pairs of energy eigenstates, and so can be treated in de- 
generate perturbation theory within the two dimensional 
subspace of the degenerate doublet. The perturbation, 
projected into the subspace of i){f) and T'tp*{r), takes 
the form: 



[^-ffsfclproj 



(33) 



where the integral runs over the unit cell and we have 
used the notation V'^(?') = ip*'^{f). Clearly, an explicit 
knowledge of the wavefunctions is needed to diagonalize 
this Hamiltonian, but a few general observations can be 
made right away. First, the eigenvalues obtained here are 
the energy splitting {SEsk) of the previously degenerate 
pair of states, which will appear in a plus-minus pair since 
the projected Hamiltonian has vanishing trace. Also, if 
we denote by Ogk the angle the vector 6k makes with the 
kx axis, the energy splitting may be written as: 



SEsk - ±VFA{0sk)\Sk\ 



(34) 



which is a massless Dirac particle with an anisotropic 
dispersion. Calculating the anisotropy function A{9sk) 
requires an explicit knowledge of the wavefunctions. 



If the vortex lattice posesses reflection symmetry about 
the y axis (x axis), then the dispersion of the Dirac nodes 
of the linearized Hamiltonian at node 1 (2) takes a par- 
ticularly simple form. The dispersion can then be written 
as: 



SE, 



'sk = ±^{v'p)Hkl + 



(35) 



where again an explicit knowledge of the wavefunctions is 
required to compute the renormalized velocities v'p and 
v'^ . The effect of the renormalization can be substantial 
and therefore the energy scale at which the Dirac node 
may be expected to make its presence felt can be quite 
different from _E2, as one might naively expect, as pointed 
out in 

Inversion Symmetry and the Dirac Node at Zero En- 
ergy: If the vortex lattice posesses inversion symmetry, 
that is, if it is invariant under the transformation r — s- — r 
(assuming the origin is the center of inversion), then it is 
easy to see from (E9|) that the superflow Pg satisfies: 



P,(-f) = -P,(f) 



(36) 



. This leads to a particle hole symmetry of the linearized 
Hamiltonian. If ip{x^y) is an eigenstate of the Hamilto- 
nian Hi ( |28| ) with energy E, then '4'{—x, —y) is also an 
eigenstate but with energy ~E [ p8[ . Inversion symmetry 
ensures there is a degenerate doublet of the linearized 
Hamiltonian Hi at the F (fc = 0) point at zero energy. 
The argument is as follows - let us focus on the spectrum 
at the F point; in which case we need to solve for the 
eigenstates of Hi on the unit cell with periodic bound- 
ary conditions - i.e. on a torus. First, consider the case 
without the Doppler term, that is a free Dirac particle 
on a torus: 

El free = VpPxCTz + VAPyCTx 

which can easily be solved. Clearly, this has a pair of 
states at zero energy, given by the product of the con- 
stant solution times any spinor. The rest of the states 
also occur in degenerate pairs, and for every pair of states 
at energy E there is a pair of states at energy —E; 
a consequence of the free Dirac Hamiltonian respecting 
the Tuirac and inversion symmetries. This spectrum is 
sketched in figure ^(a). Thus, in the free case, the spec- 
trum consists of an 'odd' number of degenerate pairs, due 
to the existence of the pair at zero energy (this can be 
made more rigorous by introducing an ultraviolet cutoff, 
and hence a finite number of states). Now, turning on the 
Doppler term Psx for inversion symmetric vortex lattices, 
preserves both the Toirac as well as particle-hole symme- 
try. Therefore we still have states coming as degenerate 
doublets, in a particle-hole symmetric spectrum. Since 
the total number of pairs of states cannot change from 
the free case, we are forced to have a degenerate doublet 
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at zero energy, so that the total number of pairs of de- 
generate states remains odd, as it was for the free case 
pof . This argument is illustrated in Figure ||. 

Notice that this argument is also valid perturbatively. 
If we imagine continuously turning up from zero the value 
of the Dopplcr term, the zero energy doublet that is 
present for the free Dirac Hamiltonian can neither split 
(it is protected by Tjjirac) nor move away from zero en- 
ergy (thanks to particle-hole symmetry) and hence con- 
tinues to exist at zero energy even for the Hamiltonian 
Hi that includes the superflow of an inversion symmetric 
vortex lattice. This can also be verified to all orders in 
perturbation theory, the details of which may be found 
in Appendix A. 

This doublet of states at zero energy for inversion sym- 
metric lattices, will give rise to a Dirac cone centered at 
zero energy, by our previous arguments. Thus, in this si- 
tation we have been able to access the nature of the low 
energy physics solely via the use of symmetry arguments. 



E 



E=0 



E 

A 

C 3 



(a) (b) 

FIG. 3. Dirac node at zero energy for inversion symmet- 
ric lattices: (a) Energy spectrum of a free Dirac particle on a 
torus (r point). There are an odd number of degenerate pairs 
of energy states, due to particle-hole symmetry and the pair 
at zero energy, (b) In the presence of the linearized Doppler 
term (Psx) for inversion symmetric lattices of ^ vortices, 
particle-hole symmetry of the spectrum is preserved, and the 
degenerate pairs of states cannot split due to symmetry under 
Toirac- To keep the count of energy levels the same as for the 
free case, there has to be a pair of states at E—0. Similar con- 
siderations apply for the case of inversion symmetric lattices 
of physical ^ vortices. 



B. Vortex Lattice of hc/2e Vortices 

We now turn to the physically relevant case of a vor- 
tex lattice of hc/2e vortices within the linearized approx- 
imation. For convenience, we collect here the relevant 
formulae derived in the previous section; the linearized 
Hamiltonian at node 1 is: 

^ vf{{p^ + a^)(T^ + {py + ayjcT^ + Psx{r)l} 

where the a fields implement the Berry phase factor for 
circling a vortex which we have taken to be: 



a = — iW log (r) 



where the are arbitrary odd integers. More generally, 
the a fields just need to satisfy: 



V X a{r) = hn'^qiS^'^^r- fi) 



(37) 



which attaches solenoids of irqi ficticious flux to the 
vortices. Clearly, physical results should not depend on 
the choice of qi. However, since we have a vortex lattice, 
it will be convenient to pick a set of qi that give rise to 
a periodic a field. Imagine choosing a set of qi, so that 
the unit cell for the problem now contains n vortices. 
Periodicity of the vector potential a requires that the 
total fiux corresponding to this vector potential vanishes 
over the unit cell, i.e.: 



E 



= 



(38) 



Where the sum is over the vortices in a unit cell. Since 
the qa are odd integers, this requires that the number of 
vortices in the unit cell, n, be an even number. Therefore, 
a vortex lattice that has one physical vortex per unit 
cell, will require at least a doubling of the unit cell when 
considered in this way. Once again, the eigenfunctions 
can be labelled by crystal momenta k that takes values in 
the appropriate Brillouin zone which is now determined, 
not just by the periodicity of the vortex lattice, but also 
by the choice of qiS. 

Invariance Under Combined Tjjirac o,rid Gauge Trans- 
formation: It is easily seen that the Hamiltonian Hf^ is 
not invariant under T^nrac] i-e. under the transformation 
-0 Tip* , the sign of the vector potential a is inverted. 
However, since the flux associated with the a field is ex- 
actly TT (times an odd integer), it is possible to arrange for 
a gauge transformation that would reverse its sign, which 
when used in combination with Toirac, would leave the 
Hamiltonian invariant. Thus, if ipj^{r) is an eigenfunction 
of H^^ with energy E and crystal momentum k, then 



^^y{r)Tr^{r) 



(39) 



is also an eigenfunction of the Hamiltonian with energy 
E. The factor of $|^| performs a gauge transformation 
that invert the sign of the gauge fields a. It is crucial that 
<i>^(r) is a single valued function, so this gauge transfor- 
mation is only allowed because we have the special case 
of TT flux solenoids. Also, it is easily shown that $^(r) is 
a periodic function, so the wavefunction 4>_j:{r) carries 

crystal momentum —k. Thus, once again, states with 
crystal momentum k and ~k have the same energy. At 
those points in the Brillouin zone that are taken to them- 
selves under time reversal, modulo a reciprocal lattice 
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vector, [k = —k (mod G)] the energy levels appear as 
degenerate doublets. This follows from the fact that the 
two wavefunctfons (V'fc(^) and a-re always 

othogonal from the antisymmetry of r. These degen- 
erate doublets will lead to Dirac cones in their vicinity 
by the same argument as before. Hence, for any vortex 
lattice we expect the spectrum of to poscss Dirac 

cones (band touchings). 

As we have noted, we are only allowed to make the 
gauge transformation above because we are dealing with 
{2m + l)7r solenoid fluxes of the ficticious gauge field 
a, that leads to ^^^^{r) being a single valued function. 
For a general value of the ficticious flux, time reversal 
symmetry of the Dirac equation is broken in an essential 
way and cannot be fixed by such a gauge transformation. 
This also implies that in general, the doubly degenerate 
states that we obtain are not accessible within perturba- 
tion theory starting with a free Dirac equation. In doing 
perturbation theory we are implicitly assuming that we 
can gradually crank up from zero the value of the ficti- 
cious flux - however since the doublet only appears at the 
TT flux values, it is missed in perturbation theory (unless 
some additional symmetry, that preserves the degenerate 
doublets for all values of the ficticious flux, is present). 
This issue is discussed in more detail in Appendix A, and 
some subtle aspects of the gauge transformation that we 
have just used are considered in Appendix B. 



a{-x,-y) = a{x,y) 



(40) 




1,0 



-1|0 



(a) 



- d - 

(b) 



FIG. 4. Different F-T gauge choices (g^) for a vortex lat- 
tice of two physical vortices per unit cell, (a) Smallest unit 
cell possible, under inversion the q and —q vortices are inter- 
changed, (b) Doubling of the unit cell due to gauge choice. 
Once again the ficticious gauge field fiux has been chosen, for 
convenience, to map onto its negative under inversion about 
the origin (*). 

Inversion Symmetry and the Dirac Cone at Zero En- 
ergy: Consider a vortex lattice that posesses inversion 
symmetry (about the origin, say) i.e. it is invariant un- 
der the transformation r — > — r. We argue here that in 
this case there exists a degenerate doublet, and hence a 
Dirac cone centered at zero energy. 

For such a lattice it is always possible to make a gauge 
choice of the set {qi\ for which, under inversion, qi — > 
—qi. A couple of examples are shown in Figure § In 
case we begin with a set of qi that do not satisfy this 
condition, one can always make a gauge transformation 
that does not alter any of the physics but brings the qis 
into this form, which happens to be convenient for the 
following analysis. In this gauge we have that the a fields 
are even under inversion, i.e. 



which is easily seen from equation (37) for the curl of a 
and the fact that the position of a g vortex is taken by a 
~q vortex under inversion. Of course we still retain the 
fact that the supcrflow is odd under inversion: 



Ps{-x,-y) = -ps{x,y) 



(41) 



These facts imply a particle-hole symmetry for the 
Hamiltonian . If tl^ix, y) is an eigenfunction of 
with energy E, then Tip*{—x, ~y) is also an eigenfunction 
but with energy -E pT[ |. 

Inversion symmetry ensures that the degenerate dou- 
blet of Hamiltonian at the k=Q (F) point is at zero 
energy. The argument is as follows - for the F point we 
need to solve for the eigenstates of Hamiltonian H['^ on 
a torus. First consider the case without the Doppler term 
and the a field, that is, a free Dirac particle on a torus. 



Hfree = VpiPxO'z 



^PyfTx} 



This has a pair of states at zero energy, given by the 
constant solution and any spinor. The rest of the states 
also occur in pairs (due to Tbirac) and for every pair 
at energy E ^ Q there is a pair of states at energy -E 
(from particle-hole symmetry). Thus in the free case, as 
sketched in Figure ^(a), there are an 'odd' number of de- 
generate doublets due to the pair at zero energy (this can 
be made more rigorous by introducing an ultraviolet cut- 
off, and hence a finite number of states). Now, turning 
on the Doppler term and the gauge field a, the energy 
levels will again appear as degenerate pairs, from sym- 
metry under the combined effect of Toirac and a gauge 
transformation, as discussed earlier. For the case of in- 
version symmetric vortex lattices, particle-hole symme- 
try is preserved as well. Therefore we have the degenerate 
doublets appearing in a particle-hole symmetric fashion. 
Since the total number of pairs of states cannot change 
from the free case, where we had an odd number of pairs, 
we are forced to have a degenerate doublet at zero en- 
ergy, as illustrated in Figure ||(b) . The doublet of states 
at zero energy will give rise to a Dirac cone centered 
at £'=0 by our previous arguments. Notice that this is 
a non-perturbative argument for the existence of a zero 
energy pair of states in inversion symmetric lattices of 
■|| vortices. In general it is not possible to access these 
states within perturbation theory since they only appear 
when the flux of the gauge field a takes on values that 
are odd multiples of tt. An alternate argument for the 
existence of this Dirac node at zero energy for inversion 
symmetric || vortex lattices, that does not invoke the 
F-T transformation, is provided in Appendix B. 

Thus, for the case of an inversion symmetric lattice of 
hc/2e vortices as well, the linearized theory predicts the 
existence of a Dirac node at zero energy. 
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IV. BEYOND THE LINEARIZED 
APPROXIMATION: MASSIVE DIRAC 
QUASIPARTICLES AND QUANTIZED 
THERMAL HALL EFFECT 

We now consider the effect of tlie subdominant cur- 
vature terms (A_ff) on the spectrum obtained from the 
hnearized equations. In view of the smallness of these 
terms compared to the hnearized Hamiltonian, their pri- 
mary effect will be to lift degeneracies that are present 
in the spectrum of the linearized problem. Therefore we 
study the effect of these terms near the Dirac cones (band 
touchings) of the linearized problem. In order to simplify 
the discussion we shall consider the the situation of a 
vortex lattice of -^(double) vortices. For the physically 
relevant case of the ||vortex lattice, the discussions runs 
on very similar lines, although IS.H in that case involves 
more terms arising from the Franz-Tesanovic gauge field 



a (26). We do not present the details of that case here. 






t = 






1 








(a) 



(b) 



FIG. 5. On including the effect of curvature terms, the 
Dirac cone obtained at zero energy for inversion symmetric 
lattices (a) now acquire a small mass gap 2mD (b). The 
chemical potential for quasiparticles, that is set by the Zee- 
man energy Ez is shown here as lying within the gap. 



A. Massive Dirac Quasiparticles 



The problem at hand then is to consider the effect of: 

Ai/ = + p2(f)]rT, + P,(f)}l -I- ^p^PyCT^ 

2m 2m pp 

(42) 

the curvature terms, on the spectrum of the linearized 
Hamiltonian at, say, node 1, 

Hi = vpiPxCTz + a^^pycr^ + Psxirjl} 



We have seen that the band structure of the linearized 
problem posesses Dirac cones at special points in the Bril- 
louin zone where there is a pair of degenerate states. This 
degeneracy is a result of the invariance of the linearized 
Hamiltonian under the symmetry operation TDirac- The 
curvature terms however do not respect this symmetry 
(under Tbii-ac, AiJ changes sign) and so are expected to 
split the degenerate doublets, giving rise to a gap in the 
dispersion (see Figure This may be analyzed within 
degenerate perturbation theory, the details of which are 
below. 

Consider a degenerate doublet (V', tiJ;*) of the lin- 
earized Hamiltonian at one of the special points in the 
Brillouin zone. We may treat the effect of AH within 
degenerate perturbation theory, since the other energy 
levels at this crystal momentum are separated by ener- 
gies of order Ei or E2, that are much larger than the 
strength of the perturbation AH. On projecting AH 
into this two dimensional subspace, it can be written as: 



[AH] 



Proj 



heB _ 
( )nAH ■ r 



(43) 



where t are the Pauli matrices that act within this two 
dimensional subspace, and there is no term proportional 
to 1, since AH changes sign under Tbirac- This defines 
for us the vector Hah- The dependence on the magnetic 
field, B, is explicitly exhibited in the prefactor. Further, 
if we neglect the last term in AH (^), that arises from 
the curvature of the gap and is smaller than the other 
terms by a factor of , then ft ah defined above is only 
a function of the anisotoropy a, and the vortex lattice 
geometry - i.e. the type of lattice and its orientation 
relative to the nodes. 

We now derive the effect of the perturbation [AH]proj 
on the dispersion around the special Brillouin zone 
points. Consider making a small excursion in crystal 
momentum k {\k\ <C 2Tr/d the typical reciprocal lat- 
tice vector) from this point. As in (p^), we obtain an 
additional term in the linearized Hamiltonian 6H^ — 
kxf^z + o;~^kyax which when projected into the two di- 
mensional space of the degenerate doublet takes the form 



[SHj:]pro] VF[kxnk,^ 



a ^kyfik^ 



(44) 



where, once again there is no term proportional to 1, 
since SHj: changes sign under Toirac- The and 
are, in general, a pair of linearly independent vectors, 
that are defined by the above equations. It may easily 
be seen that the dispersion resulting from the sum of 
these projected Hamiltonians [AH]proj + [<5-ff^]proj, is 
that of a massive Dirac particle with mass 



,heB 
mo = {—)nAH 



m 



(45) 



i.e 



i.^. the component of uah perpendicular to the plane 
defined by fik^ and fiky ■ This is most easily seen for the 
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case when n^^ and ftk^j are orthogonal Then, if n^^ 
iff'ky) is the unit vector in the direction of (^/c„) we 
can write: 

[AH + SH^]proj = v'p{kx - k^)nk^ • + v'^iky ~ k^y)nky ■ f 
+mD{nk^ xfiky) -T (46) 

. The center of the dispersion is now at (fc^, ky) — - 
■ nk^,anAH ■ "fcj and the renormahzed ve- 
locities are v'p = vplfik^ \ and v'^ ~ WA|^fc„ | and the mass 
term mo = (^)nAff • [nk^ x '^^cJ• In an appropriately 
chosen basis the above projected Hamiltonian will take 
the form: 

[AH + 5Hj:]proj = v'p{k^ - k°)Ti + v'^iky - /c°)t2 

+mDT3 (47) 

. Clearly, the resulting dispersion is: 

E{k) = ±^[v'p{k^ - fc0)]2 + [v'^iky - fc0)]2 + ^2^ (48) 

i.e. that of a massive Dirac particle with mass mp), cen- 
tered at the crystal momentum (fc°, ky) as shown in Fig- 
ure |. Of course, to calculate these quantities requires 
an explicit knowledge of the wavefunction ij:. Thus, the 
Dirac nodes obtained within the linearized approxima- 
tion acquire small gaps on inclusion of the curvature 
terms. An explicit numerical calculation of these gaps 
is presented later in this section. Below, we invoke some 
well known properties of massive Dirac particles in two 
dimensions, to derive some consequences for quasiparticle 
transport in the vortex lattice state. 

B. Quantized Thermal Transport 

The low temperature transport properties of quasipar- 
ticles depends strongly on the nature of the spectrum at 
the chemical potential. For vortex lattices that posess 
inversion symmetry, the particle-hole symmetry of the 
spectrum at each node will allow us to make some pre- 
cise statements regarding the low temperature transport 
properties - therefore in what follows we specialise to the 
case of inversion symmetric lattices. 

For inversion symmetric lattices, within the linearized 
approximation, there exists a Dirac node at the F (fc = 0) 
point centered at zero energy. On including the effect of 
the curvature terms, AH, we have seen that a gap is 
induced and the spectrum near zero energy is that of a 
massive Dirac particle, that retains particle-hole symme- 
try, and whose dispersion is centered near the F point. 
All the negative energy states are then occupied (in the 
'd' particle representation of quasiparticles that has been 
adopted) since the quasiparticle chemical potential, in 
the absence of Zeeman splitting, is at zero energy. 



This situation is topologically identical to a free Dirac 
particle in two dimensions, with a mass term. In other 
words, the spectra of the two systems can be continu- 
ously deformed into one another without closing the gap 
at the chemical potential. The massive Dirac equation in 
two dimensions has the dispersion E{p) ~ ±-\/p^ + ™^ 
and the negative energy branch is completely filled in the 
ground state. There, it is well known that the zero tem- 
perature Hall conductance of the Dirac particles is quan- 

2 

tized (Txy = 5^sgn(m) where e^i is the charge associated 
with the Dirac particle iQ. Superconductor quasi- 
particles of course do not carry a well defined electrical 
charge - however, for the case of interest the component 
of the quasiparticle spin along the applied magnetic field 
is conserved, and plays the role of eu; as we have seen 
the density of the conserved 'd' particles corresponds to 
the spin density in the direction of the field. Hence, in 
the present case of quasiparticles we expect a quantized 
spin hall conductance jl^l with ep, replaced by the 
quantum of spin Ti/2 carried by the quasiparticles. Thus, 

if we define CTq = ^^^^^ , 

^ly = \^l^g<rnD) (49) 
= (50) 

for a single node. 

It is easily seen that for inversion symmetric lattices, 
opposite nodes contibute equally to the spin hall con- 
ductivity. The Hamiltonian at node I, Hi + AH , under 
the inversion operation ■01 (^ ~^ is found to be 

identical to Hi + AH. Therefore they have the same gap 
|m£)(l)| = |m£)(i)|, and since the Hall conductance is 
unchanged by the inversion (rotation by an angle of tt), 
have identical spin Hall conductances cr^y(l) = a^y{l). 
Similarly we can relate the gaps and spin Hall conduc- 
tances between nodes 2 and 2 

If no further symmetry of the vortex lattice is assumed, 
the size of the gap and the sign of the contribution to 
a^y from the nodes 2 and 2 are not simply related to 
their values at the nodes 1 and I. Thus, as regards the 
low temperature Hall transport, two scenarios present 
themselves. The first in which the contribution from 
both pairs of nodes have the same sign, in which case 
f^^y.Tot = 4CT^y(l) = 2cr§sgn(mD ( 1) ) and the second, 
when contributions from the two pairs of nodes are op- 
posite and cancel, cr^y^pot = 20-^^(1) + 2(7^^(2) — 0, and 
the spin Hall conductance is quantized to zero. In both 
cases the quantization begins to set in at temperatures 
lower than the smallest gap. These two cases are topo- 
logically equivalent to homogenous superconductors with 
pairing wavefunctions that break time reversal symmetry, 
d.j.2_y2+idj.y in the first case and dx2_y2+is, or any other 
thermal insulator, in the second case. 

When the vortex lattice posesses an additional sym- 
metry, that of reflection about an axis that bisects the 
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nodal directions (dashed lines in Figure |l]) , we can relate 
the physics at the nodes (1 1) with that at the nodes 
(2 2). We then have |mD(l)| = |to_d(2)| and the spin 
hall conductances from all the nodes add. Thus, the first 
scenario, with a^^ j,^^ = 2crQSgn(m£)(l)) will be realized. 
For lattices that weakly break this symmetry, the ex- 
act equality of the gaps |tod(1)| and |mD(2)| no longer 
holds, but the spin hall conductance is still expected to 
remain equal to 2iTQSgn(mD(l))- A finite deformation of 
the lattice is at least required to make a transition to the 
topologically distinct state with cr^^ J,^^ = 0. 

Effect of Zeeman Splitting on Quantization: So far 
we have neglected the Zeeman splitting of the quasi- 
particles by the magnetic field, which, as we have seen 
earlier, plays the role of the chemical potential for the 
'd' particles. The magnitude of this term, Ezeeman — 
^gfisB =0.7 KTesla"^ is of the same order as the ex- 
pected mass gap in the vortex state of YBCO niD ^ 
Ef/Ep =0.5KTesla^^ and has the same field depen- 
dence. Therefore, the Zeeman term plays an impor- 
tant role in determining whether or not the exact quan- 
tization of the spin Hall conductance discussed above 
will be realised in these materials. If the quasiparticle 
gap = min[|m£)(l)|, |m£)(2)|]) exceeds the Zeeman 

term |m£)| > \Ezeeman\ then, the quasiparticle chemical 
potential lies in the gap, and quantization of the spin 
hall conductance is expected, although the temperature 
below which this will set in is now given by the differ- 
ence Imul — \Ezeeman\- If, howcver, the Zeeman split- 
tling exceeds the gap, then the chemical potential for the 
quasiparticles lies in a region where extended states are 
present, and the exact quantization will be lost. (When 
disorder is taken into account, these states may be lo- 
calised and then exact quantization will be recovered). 
The question of whether the quasiparticle gap exceeds the 
Zeeman energy depends on several details of the problem 
such as the value of the anisotropy a and the vortex lat- 
tice geometry, as well as on the value of the band mass 
and g factor of the electron for the material in consider- 
ation. 

Quasiparticle Heat Transport: While so far we have 
been discussing the spin conductance of the quasiparti- 
cles, a far more accessible transport quantity in experi- 
ments is the thermal conductance. In the superconduc- 
tor, heat is transported by the quasiparticles as well as 
the phonons. The electronic part of the thermal conduc- 
tivity is usually isolated in one of the following ways. For 
the longitudinal part of the thermal conductivity {kxx), 
the phononic contribution at low temperatures varies as 
T'^ and can typically be separated out. As for the ther- 
mal Hall conductance, that is expected to arise purely 
from the quasiparticles, since phonons are not expected 
to be skew scattered by the magnetic field . 

At low temperatures, a Widermann-Franz relation be- 
tween the themal and spin conductances of the quasipar- 
ticles is expected to hold: 



Lim 



T 



/r2 



3 {h/2y 



(51) 



Thus, for the cases discussed earlier we are led to ex- 
pect, at low temperatures, a 'quantized' thermal Hall 

conductivity, in the sense Hxy/T — + ''^{■^-^) where 
n=±2,0 for the two possible cases, and (k^x/T) 0. 

Violation of Simon-Lee Scaling In fl^ , Simon and Lee 
showed that the leading contribution to the thermal Hall 
conductivity takes the scaling form: 

(52) 



ixy^T^Fxy{T/H^^^JF^) 



Clearly, the quantized Hall conductivity derived above 
does not satisfy this scaling relation. The reason for this 
is easily seen. In deriving equatio n (p2[ ) for the vortex lat- 
tice following the arguments in ]14|], we need to assume 
that the spectrum of states in the linearized theory at 
a given crystal momentum, are well separated in energy, 
and so are only weakly mixed by the perturbing curvature 
terms (AiJ). While this assumption is true over most of 
the Brillouin zone, it breaks down at the isolated points 
where doubly degenerate states are present. It is pre- 
cisely the lifting of this degeneracy by AH that leads to 
the mass term for the Dirac quasiparticles, and the quan- 
tized thermal Hall conductivity that violates the scaling 
form above. For temperatures exceeding the induced gap 
Tod ^ Ef/Ep, Simon-Lee scaling is recovered. 



C. Numerical Evaluation of Gaps: 

In order to make the preceeding discussions more con- 
crete, we present here a numerical evaluation of the gaps 
induced by the curvature terms. These results confirm 
our assumption that their effect on the spectrum of the 
linearized Hamiltonian is indeed small. We consider the 
simple situation of a square lattice of ^ (double) vor- 
tices, oriented along the nodal directions. The mass term 
induced by AH for the Dirac cone at zero energy is calcu- 
lated (at node 1) for two values of the anisotropy param- 
eter (a=l,2). Refiection symmetry of the lattice about 
the 45^^ line allows us to relate the mass term at the other 
nodes to m]j(l). The linearised Hamiltonian to consider 
is: 



Hi = VpiPxO-z 



Pyt 



^s(r)l) 



where the explicit form of Ps{r) that we use is: 
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(53) 



(54) 



where G — ^(m, n) are the reciprocal vectors for 
a square vortex lattice with lattice parameter d, which 
from flux quantization for a — vortex lattice satisfies 
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d = y/h/eB. This form for Pg may be derived from 
equation (|2^) if wc assume that (a) Pg is divergence free 
i.e. V • Ps{r) — 0, which is vahd if the superfluid density 
can be taken as uniform and (b) the magnetic field is 
constant, which is justified if the penetration depth A is 
large compared to the magnetic length d/X <C 1. 

By our previous arguments we are assured of a pair 
of degenerate states at £^ = for this linearized Hamil- 
tonian, i.e. there exists a pair of linearly independent 
wavefunctions that satisfies Hiip{r) = 0. The pair of de- 
generate states are obtained numerically, and the effect 
of AH is evaluated in degenerate perturbation theory 
within this two dimensional subspace, to yield the Dirac 
mass rriD (^5|). Thus, the curvature terms: 



AH 



1 
2m 



P CTz 



1 
2m 



Ps{r)(Tz + —P - Ps'i- + {^PxPyO-a:} 

m py 

(55) 



generates a Dirac mass at this node 

mo — m,i + m2 + ms + {m4} 

where each term m^ results from the corresponding piece 
in AH. The term m4, which is smaller by a factor of 
A/Ep, can be neglected. The remaining terms can be 
written in the following form that explicitly displays their 
dependence on the magnetic field (B) and band mass (m). 



heB 



"^1,2,3 
XTot 



Xl,2,3 
Xl 



2m 
heB 



mn = XTot' 



2m 



where the Xi ^^'^ dimensionless quantities that depend 
only on the lattice geometry and the value of the 
anisotropy a. Therefore it is convenient to present the 
numerical results in terms of them. Physical quantities of 
interest are easily related to XTot- Its magnitude sets the 
size of the gap while its sign determines the sense of the 
quantized thermal Hall effect. Thus, the gap is given by 
2|m£)| = 2 1 XTot I and is the same at all nodes due to 
the reflection symmetry about the 45*^ line of the square 
lattice being considered. The thermal Hall response from 
all these nodes add to give the quantized thermal Hall 
conductance: 

K^y/T = 2{——sgn XTot — } 
in m 

We formulate the numerical problem in reciprocal (mo- 
mentum) space, retaining (2L+1)^ reciprocal lattice vec- 
tors. The null space of Hi {Hiip = 0) at the F point 
is obtained and the Dirac mass terms calculated using 
equations ( p3p^ , p5| ). These are shown in Figure ^ for 
two values of a, (a=l,2), where XTot as well as the sep- 
arate contribution of the different curvature terms, the 



XiS are plotted as a function of L. We see that satisfac- 
tory convergence is obtained by L=16. The magnitude 
of the Dirac mass, in both cases turns out to be roughly 
|m£,| ~ 1-2|^| (which is 0.8KTesla-\ if m is taken as 
the electron mass). Thus, for fields of a few Tesla, the 
gap is much smaller than the energy scale Ei (50KTesla2 
in YBCO), which justifies our treatment of the curvature 
terms as perturbations on the linearized Hamiltonian. In- 
cidentally, the mass term in this particular example, ex- 
ceeds the Zeeman energy (0.7KTesla~^ with g=2), and 
so the chemical potential for the quasiparticles will lie 
within the gap. 

Dirac mass term with a=1 .0 



10 12 14 



Dirac mass term with a=2.0 




FIG. 6. Numerical evaluation of the Dirac mass term for 
a square vortex lattice of hc/e vortices for two values of the 
anisotropy a = 1.0, 2.0. The magnitude of the dimensionless 
quantity XTot controls the gap size, while its sign determines 
the sign of the quantized thermal Hall effect. Here, XTot, as 
well as the separate contributions Xi arising from the differ- 
ent curvature terms are plotted as a function of L, that is a 
measure of the number of reciprocal lattice vectors (2L + 1)^, 
retained in the numerical calculation. 

The overall sign of the quantized thermal Hall conduc- 
tivity is clearly seen to be a function of a, it reverses sign 
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on going from a = 1 to a = 2 {xTot takes opposite signs 
for these two values of a). Its dependence on the sign of 
— is expected, since on inverting the sign of any of these 
(e, B, m), the quasiparticles are expected to deflect in 
the opposite direction. 

Thus, given a real material whose vortex lattice struc- 
ture and microscopic parameters are known, one can fol- 
low the recipe above to (a) calculate if a quantized Hall 
conductance is expected for the clean lattice case (i.e. 
does mo exceed Ezeeman), and if so (b) the quantization 
value (±2,0 in appropriate units), and the temperature 
scale (set by Imul — \Ezeeman\) below which quantization 
is expected. 

We have already noted the similarity of the state with 
quantized thermal Hall conductivity of ±2 to the pure 
(1^2 _y2 ± ida;j, system - both are gapped and have iden- 
tical low temperature thermal conductivities. Other au- 
thors, |l^ (2^ have considered the possibility of a mag- 
netic field inducing a 'id^jj,' component in a dx2-y2 su- 
perconductor. We believe that our work for the vortex 
lattice provides a concrete realization of the general ideas 
of [|9| 1^^. Further, we have described a proceedure 
that allows us to determine when such a state is to be 
expected, and the size of the energy gap induced. 

V. DISCUSSION 

In summary, we have considered the problem of d-wave 
quasiparticles in the mixed state within the perfect lat- 
tice approximation. Within the linearized theory we find 
a massless Dirac dispersion for the low energy quasiparti- 
cles, so long as the vortex lattice posesses inversion sym- 
metry. On going beyond the linearized approximation 
and including the smaller curvature terms, a small gap 
is found to open and the low energy excitations now be- 
have like massive Dirac particles. The size of the gap is 
proportional to the applied magnetic field, and is roughly 
of order El/Ep (« 0.5 KTesla^i in YBCO). If the sign 
of the Dirac mass term (appropriately defined) is iden- 
tical at the four nodal points, then we obtain a topo- 
logically non-trivial state - i.e. one that has gapless chi- 
ral quasiparticle modes at the edge, which in this case 
are the same as in a pure 6.^2 -y2 ± iAxy superconduc- 
tor. When the chemical potential for quasiparticles lies 
within the gap a quantized thermal Hall conductivity of 
i^xy/T = ±2 is expected at low temperatures. The 
chemical potential for these quasiparticles is controlled 
by the Zeeman energy Ez , and if this is smaller than the 
gap scale mu the quantized thermal Hall effect will be 
observed at low temperatures (T <C [mu — Ez])- We also 
described a computationally simple proceedure for eval- 
uating the numerical value of the gap and the quantized 
themal Hall conductance, which depend on the micro- 
scopic material parameters as well as the value of the 



anisotropy and the geometry of the vortex lattice. 

Although in general the quantized thermal Hall con- 
ductivity can take on any integer value (see Appendix D) , 
on the basis of the energetics for this particular situation 
we have argued that the thermal Hall coefficient is quan- 
tized to one of the discrete set of values {Kxy/T — ±2, 
in appropriate units). This is a consequence of (a) the 
independent node approximation, and (b) that at each 
of the four nodes, the only low energy feature (compared 
to Ef/Ep) is a single Dirac cone centered at zero energy 
within the linearized problem for inversion symmetric lat- 
tices. For special situations where the linearized theory 
may posess additional low energy features, for example 
as reported in |l^,|ri|,|l^ for the square lattice of vortices 
at large anisotropy (a) - other integer values of the quan- 
tized thermal Hall conductance may be realized. 

Throughout, we have treated the four nodes as being 
independent, and neglected the effect of inter- node scat- 
tering. This is expected to be important when there is 
a high degree of commensuration between the reciprocal 
vectors associated with the vortex lattice, and the mo- 
menta separating the four nodes at the Fermi surface. 
In that demonstrated in |lj] the nodal points 

acquire a gap due to the formation of a quasiparticle 
density wave. However, away from such special commen- 
suration, the independent node approximation used here 
is expected to hold. 

Our discussion in this paper has been restricted to the 
clean limit but of course disorder is always present in 
the physical system. The topological character of the 
states we have been considering renders them insensitive 
to the effects of weak disorder. Further, if the quasi- 
particle chemical potential lies outside the gap, while a 
quantized thermal Hall effect would not be expected in 
the pure case, disorder could localises the quasiparticle 
states at the chemical potential stabilizing the quantiza- 
tion. However, for larger disorder strengths, a completely 
different approach that does not rely on the Bloch nature 
of the quasiparticle states will be required. With these 
caveats in mind we turn to survey some of the relevant 
experiments in YBCO. 

Low temperature thermal-Hall measurements would 
provide the most direct test of the occurance of such 
topologically ordered states in superconductors. Cur- 
rently, such measurements on YBa2Cu307 go down to 
temperatures of about 12K in a field of 14Tesla which 
is presumably still too high to observe the quantization, 
should it exist. Indeed, although the measured K^y at a 
given temperature is found to saturate for stronger fields, 
the value at this plateau scales as T^, rather than linear 
in T as would be expected if Kxy/T was quantized. How- 
ever it is an intriguing fact that the plateau value of Kxy 
at the lowest temperature measured (12. 5K) is very close 
to what would be expected from a quantized thermal 
Hall conductance of \Kxy/T\ = 2 (in appropriate units). 
While this experiment is not conclusive with regard to 
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the low temperature state of the quasiparticles, we can 
look to specific heat measurements for a signature of the 
chemical potential lying in a gap. Such low tempera- 
ture specific heat measurements down to IK in magnetic 
fields of 14Tesla in YBa2Cu307 have been reported in [|j . 
However, they show no evidence of a gap, in fact a finite 
density of quasiparticle states that scales as the square 
root of the field was found. In principle, a quantized ther- 
mal Hall condutance could still arise in such a system, if 
the states at the chemical potential are localized - in all 
cases the quantization of the thermal Hall conductance 
will be accompanied by the vanishing of the longitudinal 
thermal conductance Kxx/T — > at temperatures below 
the (mobility) gap scale. But low temperature T < O.IK 
longitudinal thermal conductivity measurements in fields 
upto STesla in YBa2Cu306.9 (|] reveals k^x/T saturating 
to a nonzero value that rules out a quantized thermal Hall 
effect in this material down to these low temperatures. 
Thus, for this case of YBCO, either the Zeeman splitting 
causes the chemical potential to lie outside the gap re- 
gion, or else the vortex lattice in this case is so disordered 
that the perfect lattice assumption we start with requires 
serious modification. Nevertheless, given the large num- 
ber of potentially different experimental systems with a 
d-wave gap that are available, it is not unreasonable to 
expect that the quantized thermal Hall effect, realized in 
the manner described, will be observed in the future. 
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APPENDIX A: PERTURBATION THEORY 

As pointed out in [ pl| , perturbation theory may be suc- 
cessfully employed in proving the existence of zero energy 
states of the linearized Hamiltonian. Here we consider 
two cases within this technique. First, we look at general 
inversion symmetric vortex lattices of he/ e (double) vor- 
tices, where our expectation of a degenerate pair of zero 
energy states is confirmed within this technique. Sec- 
ond, we consider a specific example of a vortex lattice of 
hc/2e vortices, which is invoked later to illustrate some 
of the subtle features of the Franz- Tesanovic transforma- 
tion when used within the linearized approximation. 

Zero Energy States for Vortex Lattice ofhc/e Vortiees: 
Here, we consider an inversion symmetric vortex lattice 
of hc/e (double) vortices and establish the existence of 
a pair of zero energy states to all orders in perturbation 
theory. Essentially this is the problem of a Dirac particle 
in the presence of a periodic scalar potential V{x,y) = 
Psx{x,y): 

V{r + e,)=Vir) 

where € {61,62} are the translations that generates 
the vortex lattice and 

V{-x,-y)^~V{x,y) (56) 

which results from the inversion symmetry of the vortex 
lattice. The Hamiltonian that we consider is 

H = VF{]^ + V{x,y)} (57) 

where, by a slight abuse of notation we have found it 
convenient to define: 

To find the zero energy states we solve the equation 

H^{x,y)^0 

at the r point, i.e. with periodic boundary conditions for 
the wavefunction. We choose to work in reciprocal space, 
where, as usual, the reciprocal vectors G,-j = riigi -f- 71-252, 
with integers n = (711,712), and basis vectors (gi) that 
satisfy the conditions: 

(Ji ■ = 2Tr6ij 

The potential and the wave function may be Fourier ex- 
panded as: 



n 

T/(f)=^l7(Ga)e*^-^~ 

n 

and in this basis the eigenvalue equation for the zero 
energy states is: 

fSjiGn) + J2 ^(^« - Gn')HGR') = (58) 
ft' 

From the fact that the potential V{x, y) is odd under 
inversion (|5^), we have that the Fourier components V 
are also odd functions of the reciprocal lattice vectors: 

V{-G) = -V{G) (59) 

and in particular, the zero wave-vector component van- 
ishes y(0) = 0. 

We now prove to all orders in perturbation theory that 
there exists a pair of states that satisfies the above equa- 
tion ( psf ) . We treat V as a, perturbation on the free Dirac 
particle, and to organize the perturbation expansion it is 
convenient to imagine that there is a small control pa- 
rameter e.V^tV which is set equal to one at the end 
of the calculation. Then, if we can expand the solution 
in powers of e we obtain: 

^{G) = MG) + eMG) + e^MG) + - (60) 

. Inserting this in equation ( |5^ ) we can solve to each 
order in e to obtain, for G 7^ 0, the following recursion 
relation: 

f;4>M{G) = - ^ F(G - G')^Af-i(G') (61) 

G' 

where we have used a more compact notation replacing 
Gfi,Gfi' by G,G'. The equation for G = imposes on 
the solution the condition: 

= ^T>(-G')^M-i(G') (62) 

G' 

We start by finding t/^q, i.e. solving the free problem, 
and then use the recursion relation (^ij) to iteratively 
obtain ipM- A zero energy solution exisits if the condi- 
tion ( |6^ ) is satisfied at all orders of perturbation theory, 
i.e. if ( |62| ) is satisfied by the ipM so obtained, for all 
M. The solution to the free problem is trivial, there 
are two zero energy states corresponding to the spatially 
constant function times any spinor, namely ipo{G) = 
if G 7^ and 7Ao(0)= (1 0)^ or (0 1)^. To this order of 
perturbation theory, the condition ( |6^ ) is satisfied, since 
1^(0) = 0. Now, examinining the recursion relation ( |6l| ) 
and using (|5^), it is easily seen that at any order in per- 
turbation theory the wavefunction is an even function of 
G - that is 
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V^Af(G) =V^m(-G) (63) 

and hence the condition (|6^) for the existence of a zero 
energy state is automatically satisfied at all orders of per- 
turbation theory: 

Y,V{-G')i^M-i{G') = i^[y(-G')V'M-i(+G') 

G' G' 

= 

In fact, the two choices of ip^ give rise to two distinct 
zero energy states and -0) that are related by (t>{G) — 
Ttfj*{—G), which are easily seen to be orthogonal. 

Thus, we have been able to confirm to all orders in 
perturbation theory that the linearized Hamiltonian for 
an inversion symmetric vortex lattice of hc/e vortices 
posesses a pair of degenerate states at zero energy. 

Zero Energy States of a hc/2e Vortex Lattice: We now 
consider a vortex lattice of hc/2e vortices, specifically, 
we look at a vortex lattice with two vortices per unit cell 
but consider making a Franz- Tesanovic transformation 
that artificially doubles the unit cell, as shown in Fig- 
ure |b. The reason for considering this particular case 
is that it will help us expose some subtle features of the 
Franz- Tesanovic transformation when used in combina- 
tion with the linearized approximation, in the absence 
of a proper regularization. While we shall find a pair 
of zero energy states for this case, numerical work on 
an essentially equivalent situation ||ll| - the same lattice, 
only treated using a Franz- Tesanovic gauge transforma- 
tion that preserved the two vortex unit cell (like in figure 
^) - found a gap at zero energy. This paradoxical sit- 
uation is taken up in more detail in the following two 
appendices where a spurious coupling present in the un- 
regularised (continuum) theory is found to be responsible 
for this state of affairs. 

The linearized problem is that of a Dirac particle in the 
presence of a periodic scalar potential V{x, y) = Psx {x, y) 
and a vector potential a (chosen to be periodic) corre- 
sponding to solenoids of flux ±7rg located at the positions 
of the vortices f± ± = ± | ± | as shown in Figure ||. Thus , 
the linearized Hamiltonian we consider is: 

H = VF{{i>+^) + V{x,y)} (64) 

where the vector potential a satisfies: 

V X a(f) = TThYy^r- Rn - r++) - Sir- i?„ - f+_) 

n 

~S{r~ Rn-r-+) + 5{r- Rn+r^^)] (65) 

due to the inversion symmetry of this vortex lattice we 
have: 

= -V{r} 



and can choose the vector potential so that 

a(— r) = —a{r) 

Also, as a consequence of the artificial doubling of the 
unit cell we note the relations: 

V{r + xd)^Vir) (66) 
d{r + xd) = —d{r) (67) 

To find the zero energy states we have to solve the 
equation 

H4,{x,y)^0 

at the F point, i.e. with periodic boundary conditions 
for the wavefunction. Denoting the Fourier components 
of the various quantities by a tilde as before, we have the 
equivalent equation in reciprocal space: 

f!^{G) + Y^i^iG - G') + V{G - G')]V^(G') = (68) 

G' 

The fact that the potentials are odd functions of f 
results in: 

V{-G) = -V{G) (69) 
a(-G) = -a(G) (70) 

We consider both the vector and scalar potentials as 
perturbations of order e; V eV and a ea, and ex- 
pand the wavefunction in powers of e: 

oo 

^(G) = ^""MG) 

M=0 

. Inserting this in equation ( |5^ ) we can solve to each 
order in e which obtains for us (for G ^ 0) the following 
recursion relation: 

)G^m(G) - -^[^(G-G') + V(G-G')]^M-i(G') 

G' 

(71) 

while the equation for G = imposes on the solution the 
condition: 

= ^[^(G - G') + V{G - G')]^M-i (G') (72) 

G' 

Once again, we start the perturbation theory by find- 
ing "00, i-e. solving the free problem, and then use the 
recursion relation ( [7l| ) to iteratively obtain ipM- A zero 
energy solution exisits if the condition ( [72| ) is satisfied at 
all orders of perturbation theory. The solution to the free 
problem is trivial, there are two zero energy states cor- 
responding to the spatially constant function times any 
spinor, namely iPo{G) = if G 7^ and ipo{0)= (1 0)^ 
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or (0 l)"^. To this order of perturbation theory, the con- 
dition ( [7^ ) is satisfied, since by equation (|69|), V^(O) = 
and a(0) = 0. Now, examining the recursion relation 
( [71] ) and using (|69|), it is easily seen that at any order in 
perturbation theory the wavefunction is an even function 
of G - that is 



i^M{G)^^M{-G) 



(73) 



and hence the condition (|7^) for the existence of a zero 
energy state is automatically satisfied at all orders of per- 
turbation theory: 

^[^(-G') + y(-G'fc-i(G') = 



In fact, the two choices of ■(/'o give rise to two dis- 
tinct zero energy states {4> S'Hd ip) that are related by 
(j){r) = Ttp*{r + xd). That these wavefunctions are lin- 
early independent is shown in the next appendix. 

Thus, for this situation involving hc/2e vortices, we 
have shown to all orders of perturbation theory that 
the (continuum) linearized Hamiltonian in the particular 
Franz- Tesanovic transformation used above, has a degen- 
erate pair of zero energy states. 



APPENDIX B: REGULARIZING THE 
LINEARIZED THEORY 

In this appendix we consider some subtle features of 
the Franz- Tesanovic transformation that arise in the lin- 
earized approximation. Recall, that the ficticious U(l) 
gauge field a was introduced only to take into account 
the statistical interaction of the quasiparticles and hc/2e 
vortices. This led to the Dirac Haniiltonain: 

-fff^ = vf{{Px + ax)crz + a^^{Py + ay)crj, + x ■ Psir)} 



(74) 



with 



V X d{r) 



where qi is an odd integer. Now, as is well known from 
the theory of Dirac equations that a gauge field min- 
imally coupled to a Dirac particle has not only an orbital 
coupling, that gives rise to a phase difference for different 
paths, but in addition automatically generates a 'Zee- 
man' coupling of the 'spin' of the Dirac particle to the 
curl of the gauge field. Here, we shall refer to this cou- 
pling as the psuedo-Zeeman coupling since of course the 
physical spin of the quasiparticlc is not involved. Rather, 
the two component structure of the quasiparticlc wave- 
function arises from electron-hole mixing. 



Thus, in the linearized theory, in addition to generat- 
ing the requisite statistical interaction between the vor- 
tices and quasiparticles, we have automatically produced 
a psuedo-Zeeman interaction of the quasiparticlc with the 
curl of the ficticious vector potential a. 

Thus, the linearized Hamiltonian ( [T^ is not a faith- 
ful representation of our original problem. In addi- 
tion to the requisite statistical interaction, the Franz- 
Tesanovic transformation automatically hands us a spuri- 
ous psuedo-Zeeman interaction of the quasiparticles with 
the curl of the vector potential a, which has the form of 
delta functions of strength gjTr located at the vortex cen- 
ters. As a conseqence, this continuum theory (defined 
without an appropriate short distance regularization) is 
not invatiant under different Franz-Tesanovic transfor- 
mations, i.e. different choices of the qi. This feature has 
in fact been observed in the numerical studies of 
The origin of this problem can be traced back to the sin- 
gular gauge transformation required to turn the multiple 
valued wavefunction with branch cuts, to a single valued 
wavefunction. The transformation is singular at the vor- 
tex centers, and since the quasiparticlc amplitude (in the 
continuum theory) is finite at these points the essential 
physics of the problem is altered by this singular map- 
ping. In addition, due to these psuedo-Zeeman couplings, 
the linearized continuum theory breaks Toirac m an es- 
sential way, and would give rise to spurious gaps in the 
band structure of the quasiparticles. We believe this to 
be at the root of some numerical results such as the gaps 
reported in ||T^ for inversion symmetric non-Bravais lat- 
tices. This issue is discussed in more detail in the next 
Appendix. 

Clearly, a proper regularization of the problem that 
takes into account the vortex cores, will avoid this prob- 
lem. The simplest regularization proceedure would be 
to imagine that the quasiparticles live on a microscopic 
lattice with the vortices at the plaquette centers - then 
the quasiparticles are oblivious to the existence of the 
fiux (V X a) of the ficticious gauge field, except for the 
Aharanov-Bohm phases that they generate. 

At this stage it may appear that the conclusion ob- 
tained in Section 
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of the existence of Dirac nodes 
for the hc/2e vortex lattice case would also be suspect, 
since this was argued witin the continuum theory where 
the spurious psuedo-Zeeman coupling is expected to be 
present. However, a closer look at the argument reveals 
that we have implicitly turned off this spurious coupling. 
This is assumed when making the gauge transformation 
that combined with Toirac ( |39| ) leaves the Hamiltonian 
invariant and protects the Dirac cones. This gauge trans- 
formation is only allowed if the quasiparticles are obliv- 
ious of the psuedo-Zeeman coupling. Therefore, the re- 
sults obtained with this argument actually applies to a 
properly regularized theory (such as with a microscopic 
lattice) and reflect the physics of the original problem. 
Since the existence of the Dirac cones in the linearized 
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theory of the /ic/2e vortex lattice is one of the central 
results of this paper, we provide below an alternative 
argument that leads to this conclusion, but avoids the use 
of the Franz-Tesanovic transformation. The linearized 
Hamiltonian for a vortex lattice of hc/2e vortices before 
making the F-T transformation: 

Ti-l vfIpxCTz + a^^PyCTx + Psx{r)} (75) 

where the wavefunctions ^' are not single valued, but 
have a set of branch cuts originating from the vortex cen- 
ters, across which the wavefunction changes sign. Now, 
it is clear that Toirac is a symmetry of this problem, 
since if '0' is an eigenstate of the Hamiltonian with en- 
ergy E, then Tijj'* is also an eigenstate with the same 
energy E and has exactly the same set of branch cuts 
as i})' . Notice that this symmetry is only valid because 
the discontinuity across a branch cut corresponds to a 
sign change, and not a general phase factor that would 
be altered by complex conjugation. Further, if we select 
the branch cuts between vortices in a periodic fashion, 
we obtain a unit cell with an even number of vortices. 
Since the superflow is also a periodic function with the 
same periodicity as the vortex lattice, the eigenstates of 
T-l\ can be labelled by the a crystal momentum (fc) that 
takes values within a Brillouin zone that is determined 
not only by the vortex lattice but also by the choice of 
branch cuts. Once again, due to Toirac there is a degen- 
eracy of states with crystal momentum k and —k and for 
the special points in the Brillouin zone for which these 
two are the same, i.e. k = — fc(mod G), there is a pair of 
degenerate states. This leads to Dirac cones centered at 
these points in the Brillouin zone. For inversion symmet- 
ric lattices, it is always possible to select a set of branch 
cuts that are taken to themselves on inversion. Consider, 
for convenience, working with such a set of branch cuts 
- now it is easily seen that if i^'lr) is an eigenstate of 
the Hamiltonian Ti\ with energy i?, then ijj^—r) is also 
an eigenstate, with precisely the same set of branch cuts 
as ip'ir), but with energy —E. Thus, for inversion sym- 
metric lattices the spectrum at each node is particle-hole 
symmetric. 

We can now argue for the existence of a zero energy 
doublet at the fc = (F) point for this case of inversion 
symmetric lattices. Let us consider for definitcncss the 
case of an oblique lattice with two vortices per unit cell, 
with the branch cut chosen to go between them as shown 
in Figure (|^c) . Starting with a free Dirac particle on the 
torus (i.e. the unit cell with periodic (fc — 0) boundary 
conditions) we note that the spectrum contains a pair 
of zero energy states corresponding to the spatially con- 
stant spinors. Now, imagine introducing a hc/e (double) 
vortex at the center of the unit cell as shown in Figure 
) - as argued previously from the Tjjirac and particle- 
hole symmetry for this case, a pair of states have to still 



exist at zero energy. We can now imagine this double 
vortex being drawn apart into two hc/2e vortices with 
a branch cut running between then as shown in Figure 
(^). Since we have shown that Tjjirac and particle- hole 
symmetries continue to exist for this situation, we con- 
clude once more that a pair of zero energy states exist for 
a vortex lattice of hc/2e vortices, so long as it posesses 
inversion symmetry. This pair of states at zero energy is 
also required to be present to keep the count of the total 
number of energy states the same as in the free Dirac 
case, where there were an odd number of doublets. 




(a) (b) (c) 



FIG. 7. Alternate argument for a degenerate pair of states 
at zero energy, for inversion symmetric vortex lattices of hc/2e 
vortices. Starting with (a) a free Dirac particle on the torus 
for which a pair of zero energy states exists, it can be argued 
from the presence of Tijiracand particle-hole symmetry that 
(b) this pair persists at zero energy when a hc/e vortex is 
introduced, and (c) continues to exist at zero energy when 
this double vortex is pulled apart into a pair of hc/2e vortices, 
with a branch cut running between them. 

In this appendix wc have explored some subtle aspects 
involved with using the Franz-Tesanovic transformation 
in the linearized theory. A distinction needs to be drawn 
between a properly regularized theory, that captures the 
physics of the original problem and which is the focus of 
our work, and the naive continuum version of the theory 
that includes the psuedo-Zceman coupling and is hence 
not a faithful rendering of the original problem. As a 
consequence of the psuedo-Zeeman coupling, this con- 
tinuum theory is not invariant under different choices of 
F-T transformations and breaks T^nrac in an essential 
way and so does not in general posess Dirac cones in its 
spectrum. 

APPENDIX C: COMPARISON WITH 
NUMERICAL WORK 

In this appendix we compare the results we have ob- 
tained from symmetry considerations of the linearized 
theory, with the detailed numerical results reported in 
I p^ . From our analysis we expect to find Dirac cones 
at the special points in the Brillouin zone for which fc and 
— fc are equivalent, and further, for inversion symmetric 
lattices we expect a Dirac cone centered at zero energy. 
As we shall see these expectations are not quite borne out 
by the numerics. We suspect this to be due to the fact 
that some of the numerical work is inadvertantly study- 
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ing the naive continuum theory with spurious psuedo- 
Zeeman couphngs. In order to substantiate this claim 
we consider in more detail the naive continuum theory, 
despite the fact that it is not a faithful representation of 
the original physical problem. As we will see, a detailed 
correspondence between some of the numerical results 
and the spectrum of the naive continuum theory can be 
made. 

As has been noted, the naive continuum limit of the 
F-T transformed linearized theory involves a spurious 
psuedo-Zeeman coupling that breaks Toi^ac in an essen- 
tial way and hence, we do not expect to find Dirac cones 
(band touchings) in this theory. However, in the presence 
of other symmetries, some Dirac cones may survive. 




(a) 



(b) 
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(c) 



(d) 



FIG. 8. (a) The Brillouin Zone for a square vortex lattice, 
(b) A unit cell with two vortices as considered in reference 
[11]. (c) and (d) unit cells as considered in reference [12] 

Consider for example a square lattice of hc/2e vortices 
within a F-T gauge choice that leads to a unit cell with 
with two vortices as shown in Figure (|b) (and as con- 
sidered in ||ll[]). The Hamiltonian for this situation is: 



-fff^ = vf{{Px + ax)cTz +a ^{py+ ay)(Tx + x ■ Ps(,r)} 
with the solenoid fluxes: 

V X a{f) — Tit: q[5{r ~ dim + — )x — d{n + -^y) 

-5[f- d{m - i)x - d{n - i)y)] 



Now, while Toirac is not a symmetry of this problem, 
due to the doubling of the unit cell, there is a generalized 
'glide' symmetry - a combination of Toirac and a trans- 
lation by half a special lattice vector - that is present. 
Thus, consider the transformation. 



(76) 



. The Dirac time reversal operation reverses the sign of 
the vector-potential a, but the translation by half a diag- 
onal lattice vector interchanges the q and —q fluxes and 
hence restores invariance. Of course, the superflow is un- 
affected by this translation as it is periodic with a smaller 
unit cell. This leads to a degeneracy between energy lev- 
els with crystal momentum k and — fc and for those points 
in the Brillouin zone that are taken to themselves under 
this operation (labelled F, X, Y and M in Figure ||a). At 
these points we have a pair of wavefunctions tpki^ ^^'^ 
T7/>l(f + ^{x + y)) with the same energy eigenvalue. In or- 
der that we have a degenerate doublet, we need that this 
pair be linearly independent. Earlier, when the symme- 
try operation was simply Toirac, the linear independence 
followed from the fact that the two wavefunctions were 
othogonal from the antisymmetry of t. Here however we 
will need to adopt a more indirect line of reasoning, par- 
alleling that used to establish Kramer's degeneracy for 
electrons |33 . To learn whether these wavefunctions are 
indeed linearly independent at these special points in the 
Brillouin zone (for k e {0, ^x, ^y, ^{x + y)} the F, 
X, Y and M points), we adopt a proof by contradiction. 
Assume that the pair of functions are linearly dependent. 
Then, since they have the same normalization, they can 
at best differ by a phase factor e*^. 



^fc(r1=e^''TVI(r+-(x + y)) 



(77) 



Performing the 'glide' symmetry operation on both 
sides of the equation we obtain: 



^^j:{r + d{x + y)) (78) 



where we have used the fact that = —1. Using now 
the fact that translation by d{x + y) amounts to a lattice 
translation within this expanded two vortex unit cell 

i^^{r + d{x + y)) = e'~^<-d+yd)^^^^^ 

and substituting equation ( |77| ) for the left hand side of 
([78|) we obtain: 



(79) 



which for the X and Y points [k S {^x, ^y}) yields the 
consistent equation ipj^ir) = '0j('^) whereas for the F and 
M points (fc e {0, -I- fj/}) this yields a contradicition 
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tp^{r) = —tp^{r). Thus we can conclude that the pair of 
states at the T and M points have to be hnearly indepen- 
dent and the spectrum at these points in the Brillouin 
zone consists of degenerate doublets. For the X and Y 
points, this conclusion cannot be drawn. As a conse- 
quence, for the situation in Figure ||b, we expect Dirac 
cones at the F and M points but not necessarily at the 
X and Y points, in the continuum theory band struc- 
ture. Further, inversion symmetry of this vortex lattice 
gives us a particle hole symmetric band structure - if 
ijjir) is an eigenstate with energy E, then T'ip*{—r) is an 
eigenstate with energy —E. As we have reasoned several 
times before, the particle hole symmetry combined with 
the doubly degenerate nature of the spectrum at the F 
point leads to a Dirac cone centered at zero energy. This 
is exactly the result of the numerical studies in 1 11 1 for 
this case. Dirac cones are found at the F and M points, 
but not at the X and Y points of the Brillouin zone and at 
the F point there is a Dirac cone centered at zero energy. 

A similar analysis for the situation of Figure (||c and d) 
can be carried out, and the continuum theory in these two 
cases are expected to have Dirac cones at the F, X and 
M points and F, X and Y points respectively as well as 
a Dirac cone centered at zero energy in both cases. This 
is exactly as seen in the numerical work on the linearized 
model in 1131 . 
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(a) (b) 

FIG. 9. An oblique vortex lattice of hc/2e vortices with two 
vortices per unit cell. In (a) a F-T transformation that pre- 
serves this periodicity while in (b) a F-T transformation that 
doubles the unit cell. Case (b) can be shown to have a pair of 
zero energy states even in the continuum (unregularized) the- 
ory, that are protected by a 'glide' symmetry resulting from 
the artificial doubling of the unit cell. 

If instead we considered a case where such a 'glide' 
symmetry was absent, then in the continuum theory we 
would not obtain Dirac cones anywhere in the Brillouin 
zone. For example if we have a vortex lattice that con- 
tains two vortices per unit cell, and make an F-T trans- 
formation that preserves this periodicity, as for example 
in Figure ^ we do not expect to find Dirac cones and 
there will be no states at zero energy. This is also the re- 
sult obtained in for the case of a non-Bravais lattice 
with two vortices per unit cell, within a F-T transforma- 
tion that preserves this periodicity. If however we were 
to consider this same problem but within a F-T trans- 
formation that artificially doubles the unit cell, as shown 
for example in Figure then a 'glide' symmetry will be 
obtained. By the arguments in this Appendix we expect 
to find Dirac nodes and a pair of zero energy states in 



this continuum theory. The existence of the zero energy 
states in this case was explicitly verified in Appendix A 
within perturbation theory, where the 'glide' symmetry 
protects the zero energy states for all values of the flux 
q (not just odd integers) and allows it to be accessed 
within perturbation theory. This is a clear indication of 
the non-invariance of the continuum theory to different 
Franz- Tesanovic gauge transformations. 

Such non-invariance of the spectrum to different FT 
transformations was reported in ]l^ where a reciprocal 
space approach was used to numerically diagonalize the 
linearized Hamiltonian. Clearly, in this scheme it is the 
continuum theory which is being studied. Thus, the un- 
physical psuedo-Zeeman coupling will play a role and fur- 
ther, due to the cutoff in momentum space the solenoids 
of ficticious flux will acquire a finite thickness and have 
an effect beyond just implementing the Aharanov-Bohm 
phase factor of (-1). Hence the non-invariance with re- 
spect to different F-T transformations and the detailed 
agreement with our expectations for the continuum the- 
ory. On the other hand the tight binding lattice approach 
also reported in the same work , which does not adopt 
a linearization was, naturally, found invariant to different 
choices of F-T transformations. It is less clear, however, 
as to why the real space diagonalization approach of |Tl| ] 
also appears to be simulating the continuum theory. We 
suspect this could arise if the solenoids of flux happenned 
to be smoothed out over over a few plaquettes - which 
would have the same qualitative effect as the presence of 
a psuedo-Zeeman coupling. In principle however, such a 
real space approach should be able to provide an appro- 
priate regularization for the F-T transformed linearized 
problem. One particularly convenient way to implement 
the (-1) phase factors would be to use Ising link variables. 

In summary, for the F-T transformed linearized theory, 
it is important to draw a distinction between a properly 
regularised theory and the naive continuum limit. The 
latter is plagued by the presence of the psuedo-Zeeman 
coupling that renders it non-invariant under different F- 
T transformations, breaks Toirac in an essential way and 
is not a faithful representation of the original problem. 
We suspect that some of the results reported in the pio- 
neering numerical work on this linearized approximation 
|]lT| fl^ may have actually been studying this continuum 
theory. In order to substantiate this, in this appendix we 
have analysed the continuum theory using the methods 
of this paper, and drawn conclusions regarding the na- 
ture of the spectrum, that are found to agree well with 
some of these numerical results. 



APPENDIX D: CHERN NUMBERS FOR 
QUASIPARTICLE BANDS 

Consider calculating the spin Hall conductivity of 
the full quasiparticle Hamiltonian Hsdo (the linearized 
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Hamiltonian has vanishing spin Hall conductance due 
to invariance under T^irac)- We use the Kubo formula 
that evaluates the response of the quasiparticle current 
jj, = '"^ to a perturbation corresponding to a uniform 
Zeeman gradient (e) SH = eyl along the orthogonal di- 
rection. 

At zero temperature the spin Hall conductivity for an 
isolated filled band takes a particularly simple form, just 
as in the case of the electrical Hall effect with charged 
particles in a periodic potential [ p4| . It is proportional 
to a topological invariant, the Chern number |9j that can 
be asigned to the quasiparticle band labelled by 'n' via 
the formula: 

<"N = ^(|)'C^n (80) 

Cn = ^ f Vfe X A{k)d^k (81) 

The vector potential is defined in terms of the wavefunc- 
tions of the quasiparticles in that band by 

f (Pr 

A{k)=-i j — {<fc(?')VfcM„fc(0 +<fcWVfeW„fc(r)} 

(82) 

where the integral is over the unit cell, and we have 
written out the components of the quasiparticle wave- 
functions of band 'n' as ipnkir) — (u„fc(r) u„fc(r))^ that 
we choose to be normalized according to: 

1=/ ^(Kfc(r-)P + Kfc(r)p) (83) 

J unit — cell ^ 

where V is the area of the unit cell. Clearly, a phase redef- 
inition of the wavefunction at each point of the Brillouin 
zone V'nfc e^'^^'^^'^nk corresponds to a gauge transfor- 
mation A{k) A{k) + Vfex(^)- This does not affect 
the value of the Chern number - which is thus defined 
despite the phase ambiguity of the wavefunctions. Fol- 
lowing standard arguments developed for electrons [ p5| , 
it is easily established that for the case of superconduc- 
tor quasiparticles as well the C„ are integers and hence 
the spin Hall conductivity, in appropriate units, will be 
quantized to integers for this case of quasiparticles of a 
spin-singlet superconductors, so long as the chemical po- 
tential for quasiparticles lies in a gap For the 
case of homogenous superconductors that break Time re- 
versal and Parity symmetries, the appropriate limit of 
experssion (^,^2|) reduces to the topological invariant of 

The quasiparticles therefore, could be said to blend 
together Bloch wave and Landau level properties. While 
they form bands that are labelled by Bloch crystal mo- 
menta, these bands, like Landau levels, carry a Chern 
number that is related to the spin-Hall conductivity. 
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